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Abstract 



' We consider the role of hadronic rescattering of daughter kaons on the 

in 

observed mass spectra from (f> meson decays in ultra-relativistic heavy ion 

■ collisions. A hadronic cascade code (RQMD v2.4) shows that ~26% of all 

ON " 

0's decaying to K + K in central Pb+Pb collisions at SPS energies {Ebeam = 
158GeV/A) have a rescattered or absorbed daughter. This significantly affects 
the reconstructed invariant mass of the pair and shifts (j) mesons out of the 

a > 

mass peak. Kaon rescattering depletes the low velocity region, hardening 
^ ■ and broadening the observed phi mt and rapidity distributions respectively, 



relative to the dilepton channel. This effect produces an apparent change 
in the experimentally determined branching ratio not necessarily related to 
chiral symmetry restoration. Comparisons to recent experimental measures 
at CERN energies reveal a possible mechanism to account for the shape of 
the observed spectra, though not their absolute relative magnitude. 
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A. Introduction 



Chiral symmetry, the symmetry which acts separately on left- and right-handed fields 
in QCD, is spontaneously broken in nature, resulting in associated Goldstone bosons in the 
pion field as well as the observed hadron mass spectra. However, it is expected that the 
spontaneously broken part of chiral symmetry may be restored at high temperatures and 
densities |fl. Such conditions might be reached in the interior of neutron stars or during a 
relativistic heavy ion collision. In fact, a primary goal of high energy nuclear collisions is 
the creation of such matter and the study of its properties. 

One expected signal of the restoration of chiral symmetry is a change of the vector meson 
properties 0. In the nuclear system, the self energy of hadrons is changed by the medium 
they inhabit. A change in the effective mass of daughter particles can change the effective 
lifetime, and consequently the observed width, of the parent particle that decays in medium. 
Such modifications could be accompanied by a measureable change in the branching ratio 
of mesons decaying in the medium. The light vector mesons (p, to, 0) offer an especially 
promising channel for study of chiral symmetry effects due to their multi-channel decays and 
lifetimes comparable to the space-time extent of the system produced in heavy ion collisions. 
The (f) is of particular interest as the sum of the daughter mass of its di-kaon decay channel 
is very close to the mass of the 0. As a result, even a small change in daughter or parent 
masses could measureably alter the decay channel 0. 

Study of vector mesons through their leptonic decay channels, either e + e~ or 
should provide the cleanest signal of the changing masses as leptons interact with the nuclear 
medium predominantly electro magnetically. In contrast, decays to hadrons are affected by 
strong final state interactions. If a decays in the center of the reaction volume in a heavy 
ion collision, the lifetime of the kaon travelling a distance d through a length L of hadronic 
matter is exp(-Lda). Approximating the heavy ion collision as a pion bath with density 
d = .5/fm 3 and cross section a = 10 — lOOmb, the 1/e pathlength of a kaons would be 
between .2fm and 2fm, much smaller than the size of the collision region. However, while 
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the hadronic daughters interact in the medium, leptons should escape unscathed. 

The rescattering of hadronic daughters in the nuclear medium could mimic or obscure 
effects of chiral symmetry restoration by causing the reconstructed invariant mass to fall 
outside the vector meson peak, effectively decreasing the measured yields. While studies of 
the viability of the di-kaon channel of the <fi for studying chiral symmetry restoration have 
been done |2|||, the effect of the rescattering of daughters on the experimentally measured 
branching ratio has not previously been considered. 

B. Model Results 

To study the effect of hadronic scattering of daughters, we implemented the hadronic 
cascade code RQMD version 2.4 to describe the space time distribution of 0's and their 
daughters. In RQMD one can follow the history of all 0's that decay throughout the collision 
along with their daughter kaons. Upon simulation of an event we determine the positions and 
momenta of all kaons originating from decays. Fig. [I] shows the invariant mass distribution 
of all kaon pairs from <fi decays in simulated central Pb+Pb events at 158 GeV-A/c beam 
energy. The right hand figures show 0's whose daughters escaped from the collision zone 
without rescattering (top figure) and those whose daughters did rescatter (bottom) yet still 
escaped the collision zone as kaons. The left hand figure gives an overlay for comparison. 
Since the (f) lifetime is comparable to, though larger than, the expected lifetime of a heavy 
ion collision, most 0's do not decay in medium resulting in a tight peak of the invariant mass 
(Mi nv ) with a width of ~4 MeV. However, ~ 17% of 0's decaying to two kaons have at least 
one daughter kaon that rescatters and another (non-orthogonal) ~ 17% of 0's have at least 
one daughter that is absorbed. The sum of these two processes leads to ~ 26% of decaying 
0's with an absorbed or rescattered daughter. The Mi nv distribution of these kaons is much 
broader than the original distribution; broad enough to escape detection, in particular in 
experiments with low statistics or large combinatorial backgrounds. 

Tables [TJ and |TI| itemize the channels contributing to the first hadronic scattering of a 
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daughter kaon from a phi meson. The dominant interactions for both K + and K~ daughters 
proceed through the K* resonance with a pion or through annihilation of the strange quark 
by interaction with another kaon. Nearly one-quarter of all rescatterings occur through a 
variety of high mass resonances and a sizeable fraction of the rescattered K~ daughters 
proceed by a reaction with a nucleon into a X or other high mass strange baryon. 

The circles in Fig. ||] display the <fi survival probability, i.e. the probability that neither 
of the daughter kaons rescatter, as a function of m t for the rapidity region \y\ < 1. The 
probability decreases with m t , approaching ~ 60% at the lowest m t . As a result, the 
measured yield of 0's through the kaon channel should be lower at low m t than that measured 
in the dilepton channels. The ratio of the yield of <p — > K + K~ to — > corrected by the 
branching ratio should then approximately follow the circle symbols in Fig. ^| if daughter 
rescattering is the sole mechanism contributing to the difference. 

Figs. H and |4] show reconstructed yields from RQMD as a function of m t and space-time, 
respectively for \y\ < 1. Fig. |^ shows d 2 N^/mtdmtdy at mid rapidity for all <j) mesons which 
decay in RQMD overlayed with the m t distribution for those 0's whose daughters kaons do 
not rescatter and those 0's whose daughters do rescatter. The depletion of reconstructed 0's 
at low m t results in a higher effective temperature of the meson at low m t . The inverse 
slope for < m t — < 1 GeV of the original 0's is T = 220 ± 3 MeV while for 0's observed 
through the kaon decay channel RQMD predicts T = 242 ±4 MeV. The extracted yields and 
temperatures of 0's should therefore be measureably different due to daughter rescattering. 

The m t dependence of the rescattering effect upon 0's is a reflection of the phase space 
freeze-out distribution of particles in heavy ion collisions. Fig. f| shows the freeze-out distri- 
bution in time and space of the 0's from RQMD. The left hand figure shows the freeze out 
longitudinal proper time r = \Jt 2 — z 2 where t and z are the freeze-out time and z-position, 
while the right hand figure displays the radial freeze-out distribution. The solid line corre- 
sponds to those 0's whose daughters escaped unscathed from the collision, while the dashed 
line shows the distribution for those whose daughters did rescatter. It is clear from the 
figure, and intuitive, that those 0's with rescattered daughters predominantly decayed at 
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early times in the dense nuclear media, during that part of the collision of particular interest 
for chiral symmetry measurements. 

This picture is echoed in the rapidity distributions in Fig. [5] where in the top plot we 
present the rapidity of (1) all <p's, (2) those 0's whose daughter kaons escape the collision zone 
unperturbed and (3) those 0's whose daughters do not escape. Here we see the approximately 
25% of all 0's that are lost in the collision. The effect of rescattering is similar to that seen in 
the rrit distribution as those 0's closer to y = have a larger probability of being rescattered. 
In the bottom half of Fig. [| we plot the probability of survival of both daughter kaons which 
has a clear rapidity dependence. This dependence leads to a marginal widening of the 
measured — > KK rapidity distribution relative to that of the <f> — > 11 channel. 

The conclusions from RQMD simulations are intuitive. Those 0's that decay at early 
times are more apt to have rescattered daughter kaons due to the amount of hadronic 
matter through which they must travel. Further, rescattering and pressure build-up during 
the collision implies that particles that freeze-out early in the collision will have a softer 
transverse mass distribution than those which freeze-out later ||. We then expect that the 
rescattering effects at early times will be reflected in the transverse mass distribution as 
confirmed in RQMD. 

In order to compare these observations with experimental measurements we note that 
RQMD only includes the imaginary part of the cross section in its calculations. While 
Shuryak and Thorrsen [[H]] have shown that the real cross section for kaons in a dense 



nuclear medium is much smaller than the imaginary cross section, we estimate here the 
maximal effect on the observed cross section expected from the inclusion of the real part. 

The addition of the real cross section can, at most, rescatter or absorb all kaons from 
0's that decay within the freeze-out volume. Shown in the bottom two plots of Fig. |] is 
the point of last rescattering for kaons from RQMD. We define the freeze-out volume for 
kaons to be that point within which 95% of all kaons have had their last interaction, which 
corresponds to r < 36.5fm and r < 20fm. The maximal effect then of adding the real 
cross section to the RQMD calculation would be having all 0's that decay within these (r,r) 



bounds be unreconstructable and lost from the invariant multiplicity. 

The square symbols in Fig. |2] show the result of making such a drastic assumption. In 
the lowest m t bin approximately 25% more 0's are depleted than in the imaginary only 
calculation. Note that the inclusion of the real part of the cross section will only change the 
quantitative result from RQMD slightly while the qualitative shape remains the same. 

C. Discussion 

Recently, two experiments at the SPS studying Pb+Pb collisions with a beam energy of 
158GeV/c have made preliminary reports of measurements ^-0]. Experiment NA50 || 
measured ~ * at mid rapidity over the transverse mass (mt) range 1.7 < m t < 3.2 

GeV/c 2 while experiment NA49 |||7j reported a — > K + K~ distribution also at midrapidity 
but for 1 < m t < 2.2 GeV/c 2 . The reported m t inverse slopes are strikingly different; NA50 
quotes T = 218 ± 10 MeV while NA49 finds T = 305 ± 15 MeV. 

Although, within the present accuracy, the data seem to be more consistent in the m t 
range where they overlap, the extrapolated yields are significantly different. This points 
either towards a drastic softening of the m t distribution with increasing m t or distinctly 
different spectra reconstructed from — > KK and — > fifi. The latter is in qualitative 
agreement with the effect of rescattering of the decay kaons, depleting the low m t region. 
Quantitatively RQMD predicts a 17 MeV difference of the slope, much smaller than the 
observed difference. The curve in Fig. |2| shows the ratio of the two extrapolated spectra, 
— > KK/(f) — > fifi, compared to the calculations from RQMD for the expected and maximal 
effect of rescattered daughter kaons. The curve is well below what could be described by 
even the maximal daughter rescattering and we conclude that this effect can not by itself 
describe the experimental data. 

It is interesting to note, however, that RQMD does reproduce the observed — > KK 
rapidity distribution. The line in Fig. [| corresponds to the experimentalists gaussian fit to 
their data in |7J with no renormalization on our part. This line corresponds quite nicely to the 
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RQMD curve for measured <p KK for those phi's whose daughters did not rescatter. The 
gaussian width and height of the distribution from RQMD (er = .96 ± .1, A — 2.45 ± .46) are 
approximately consistent with the fit to the experimental data (cr = 1.22±.17, A = 2.43±.15) 
though the RQMD distribution is not particularly well described by a gaussian. 

Comparisons of — > KK with <f> — > 11 may be very informative if the spectra are 
measured over the same range of m t , with similar systematics. If the ratio of the m t spectra 
has the shape characterstic of rescattering, the low m t dip in the — > KK/ll ratio reflects 
the amount of rescattering and therefore the time spent in the dense hadronic phase. This 
could help to clear up uncertainties about how long the hadronic system interacts before 



freezing out |TT],[T2|. Including effects of chiral symmetry restoration on kaon properties 
may alter these arguments, but it is unlikely that both effects will produce identical m t 
dependencies. 
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TABLES 



Interaction 


Percentage 


K+ir -► K* 


24% 


K + K~ — > 7T7r or or / 


26% 


K+i^o _^ ^ or ^+A° 


20% 


A + High Mass Resonance — ► X 


27% 


All others 


3% 



TABLE I. The dominant channels in the rescattering of K + daughters from <p decays in Pb+Pb 



collisions at SPS energies. 



Interaction 


Percentage 


K~K + — > 7T7T Or 7T/9 Or TTT] 


23% 


K-K° -»• vrr/ or A A 


10% 


A^vr -► if* or Atf 


21% 


A~n^ A 


8% 


A~p — > A 


17% 


A - High Mass Resonance — > A 


21% 



TABLE II. The dominant channels in the rescattering of A daughters from <fi decays in 



Pb+Pb collisions at SPS energies. 
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FIG. 1. Invariant mass distributions of reconstructed 4> — > K + K~ in RQMD. The top right 
hand figure shows the reconstructed peak for 0's with non-rescattered daughter kaons. The bottom 
right figure is for those </>'s who had either one or both daughters rescatter before leaving the collision 
zone. The left hand figure shows an overlay of these two plots. 
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FIG. 2. The probability that a <j) produced in RQMD that decays to two kaons will have 
daughters that escape the collision zone without rescattering as a function of m t (circles) compared 
to experimental measurements of <fi — > \x\i (NA50) to <fi — > KK (NA49) corrected for the appropriate 
branching ratio [5-7] (dotted line). The squares represent the maximal possible depletion of <p 
mesons as described in the text. All points corresond to those 0's with rapidity \y\ < 1. 
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FIG. 3. The transverse multiplicity distributions of all 0's that decay to two kaons in RQMD 
(crosses), as well as those those whose daughters rescatter (squares) and those whose daughters 
do not rescatter (triangles). All multiplicity distributions have been corrected for the 4> branching 
ratio to K+K~. 
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FIG. 4. The two top plots show the space time freeze out distributions for $>'s whose daughters 
rescattered (dashed) as well as those whose daughters escaped the collision zone without rescat- 
tering (solid). The bottom two plots show the freeze out distribution of all kaons from Pb+Pb 
collisions at the SPS. r is the longitudinal proper time r 2 = t 2 — z 2 and r is the radial position 
r 2 = x 2 + y 2 . 
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FIG. 5. Rapidity distributions of 0's decaying to K + K~ corrected for the branching ratio 
from Pb+Pb collisions at the SPS determined by RQMD v2.4. The top plot shows the rapidity 
distribution for all </>'s (top), all </>'s whose daughters were not rescattered or absorbed (middle), 
and all </>'s whose daughter were rescattered or absorbed (bottom). The curve corresponds to a fit 
to the experimental data as described in the text. The bottom plot displays the survival probability 
as a function of rapidity 
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